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ABSTRACT 

The M31 globular cluster X-ray binary XB158 (a.k.a. Bo 158) exhibits intensity dips on a 2.78 
hr period in some observations, but not others. The short period suggests a low mass ratio, and an 
asymmetric, precessing disk due to additional tidal torques from the donor star since the disk crosses 
the 3:1 resonance. Previous theoretical 3D smoothed particle hydrodynamical modeling suggested a 
super-orbital disk precession period 29±1 times the orbital period, i.e. —81±3 hr. We conducted a 
Swift monitoring campaign of 30 observations over —1 month in order to search for evidence of such a 
super-orbital period. Fitting the 0.3-10 keV Swift XRT luminosity lightcurve with a sinusoid yielded 
a period of 5.65±0.05 days, and a >5cr improvement in x 2 over the best fit constant intensity model. 

A Lomb-Scargle periodogram revealed that periods 5.4-5.8 days were detected at a >3er level, with a 
peak at 5.6 days. We consider this strong evidence for a 5.65 day super-orbital period, —70% longer 
than the predicted period. The 0.3-10 keV luminosity varied by a factor —5, consistent with variations 
seen in long-term monitoring from Chandra. We conclude that other X-ray binaries exhibiting similar 
long-term behaviour are likely to also be X-ray binaries with low mass ratios and super-orbital periods. 
Subject headings: x-rays: general — x-rays: binaries 


1. INTRODUCTION 

The M31 globular cluster (GC) B158 (a.k.a. Bo 158), 
named following the Revised Bolog na Catalogue v3.4 
(jCalleti et al .1120041 12006112007112009T) . contains a bright 
X-ray s ource that was dis cover ed by the Einstein obser¬ 
vatory l|Trinchieri fc FabbianolfT991l) . and has shown up 
in every X-ray observation of the region since; we call 
this X-ray source XB158. 

XB158 exhibited strong intensity modulation on a 
10017±50 s (—2.78 hr) perio d during the 2002 Januar y 
XMM-Newton observation (ITrudolvubov et all l2002f) . 
iTrudolvubov et ahl (120021) found similar variation in the 
folded lightcurves from a 1991 June ROSAT observation 
and a 2000 June XMM-Newton observation; they found 
that the amplitude of modulation decreased with increas¬ 
ing source intensity. They found no such modulation 
in the 2001 June XMM-Newton observation, setting a 
2 a upper limit of 10% modul ation. Assuming that this 
represents the orbital period, ITrudolvubov et ahl (|2002f ) 
found that this is probably a neutron star binary with a 
low mass donor with a separation <10 11 cm (i.e. a low 
mass XB-ray binary, LMXB). 

However, analysis of the unfolded 2000 June XMM- 
Newton lightcurve revealed a single deep dip at the end 
of the observation, with n o evidence for dips in the 
two previ ous orbital cycl e s (IBarnard et al.l 1200(1 ). Fur¬ 
thermore, iBarnard et al.l (|2006t ) analyzed 3 proprietary 
XMM-Newton observations over 2004 July 17-19, finding 
—100% dipping for one orbital cycle, and zero evidence 
for dips in other cycles; we concluded that the disk was 
precessing, with dips only visible for some part of the 
super-orbital cycle. 

Such behavior is associated with the “superhump” phe¬ 
nomenon that is observed in accre ting binaries where 
the mass ratio is smaller than —0.3 (IWhitehurst fe Kingl 


119911) . Super humps were first identified in the super¬ 
outbursts of the SU UMa subclass of cataclysmic vari¬ 
ables (accreting white dwarf binaries), manifesting as a 
periodic increase in the optical brightness o n a period 
that is slightly longer than the orbital period dVoutll 1 9741 : 
!Warneilfl975f ). SU UMas are a subclass of dwarf novae 
with orbital periods <2 hr, that exhibit occasional su¬ 
peroutbur sts that las t >5 times as long as the normal 
outburst s (1 Vogtl1 1980h . 

lOsakil (|1989t ) proposed that these superoutbursts are 
enhanced by a tidal instability that occurs when the 
outer disk crosses the 3:1 resonance with the secondary; 
the additional tidal torque causes the disk to elongate 
and precess, and also greatly enhances the loss of angu¬ 
lar momentum (and therefore the accretion rate). The 
disk precession is prograde in the rest frame, and the sec¬ 
ondary repeats its motion with respect to the disk on the 
beat period between the orbital period and the preces¬ 
sion period, a few percent longer than the orbital period. 
The secondary modulates the disk’s viscous dissipation 
on this period, giving rise to the maxima in the opti¬ 
cal lightcurve known as superhumps. Some short period, 
persistently bright CVs exhibit perman ent superhumps 
(jPattersonlll999ilRetter fc Navloill2000D . 

4U 1916—053 is a high inclination neutron star LMXB 
with an X-ray period of 5 0.00U0.08 min and an optical 
period 50.458±0.003 min (|Callanan et al.lll995D . It ex¬ 
hibits periodic X-ray intensity dips due to absorption by 
material in the outer disk; the amplitude of these dips 
varies over —0 to —100% on a —4 day period , the preces¬ 
sion pe riod of the disk dChurch et al.l [19971 : iChou et all 
120011) . lHaswell et all (120011) showed that NS LMXBs 
with orbital periods shorter than —4.2 hr are likely to 
exhibit superhumps, and identified 4U 1916-053 as a per¬ 
sistent superhumping source. XB1 58 appears to be som e- 
what analogous to 4U 1916—053 (IBarnard et al.ll200 (>). 
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Table 1 

Journal of Swift observations. For each observation we give the time (in days) after the first observation, pointing position for the 
observation, XRT exposure time, and off-axis angle for XB158. We also give the net source counts from the observation 


Obs 

T — To 


RA 

Dec 

XRT EXP 

T>/' 

Counts 

00032702001 

0.0 

00 

43 

00.81 

+41 

17 

00.2 

2234 

10.0 

57 

00032702002 

0.9 

00 

43 

06.53 

+41 

16 

54.6 

2510 

9.7 

41 

00032702003 

2.0 

00 

43 

03.73 

+41 

16 

40.8 

2299 

9.5 

13 

00032702004 

3.0 

00 

43 

06.35 

+41 

16 

51.3 

2299 

9.6 

40 

00032702005 

3.1 

00 

43 

13.10 

+41 

16 

30.4 

2489 

9.2 

43 

00032702006 

4.4 

00 

43 

03.22 

+41 

13 

58.5 

2572 

7.0 

65 

00032702007 

5.4 

00 

43 

10.44 

+41 

15 

10.9 

2511 

7.9 

51 

00032702008 

6.4 

00 

43 

06.44 

+41 

16 

06.4 

2635 

8.9 

57 

00032702009 

7.3 

00 

43 

11.97 

+41 

16 

27.7 

2476 

9.1 

29 

00032702010 

8.8 

00 

43 

03.76 

+41 

16 

15.5 

2644 

9.1 

39 

00032702011 

9.5 

00 

43 

08.40 

+41 

17 

27.6 

2405 

10.2 

53 

00032702012 

10.1 

00 

43 

07.20 

+41 

16 

29.8 

2502 

9.2 

43 

00032702013 

11.2 

00 

43 

03.63 

+41 

17 

36.2 

2514 

10.5 

48 

00032702014 

12.4 

00 

43 

04.82 

+41 

17 

00.3 

2411 

9.8 

38 

00032702015 

13.2 

00 

43 

07.68 

+41 

15 

24.0 

2490 

8.2 

11 

00032702016 

14.2 

00 

43 

01.12 

+41 

15 

11.3 

2511 

8.2 

27 

00032702017 

15.1 

00 

43 

06.35 

+41 

14 

46.6 

2273 

7.6 

54 

00032702018 

16.8 

00 

43 

02.79 

+41 

17 

00.9 

2314 

9.9 

55 

00032702019 

17.1 

00 

43 

08.59 

+41 

15 

26.4 

1992 

8.2 

45 

00032702020 

18.1 

00 

43 

04.88 

+41 

17 

31.5 

1389 

10.3 

13 

00032702021 

19.4 

00 

43 

03.24 

+41 

17 

28.0 

2690 

10.3 

18 

00032702022 

20.2 

00 

43 

02.68 

+41 

17 

26.6 

2686 

10.3 

42 

00032702023 

21.1 

00 

43 

03.07 

+41 

17 

07.0 

2821 

10.0 

65 

00032702024 

22.3 

00 

43 

04.56 

+41 

16 

19.5 

2659 

9.2 

50 

00032702025 

23.2 

00 

43 

10.49 

+41 

16 

54.1 

2711 

9.6 

46 

00032702026 

24.3 

00 

43 

07.86 

+41 

15 

05.1 

2740 

7.8 

35 

00032702027 

25.1 

00 

43 

06.42 

+41 

16 

41.0 

2424 

9.5 

18 

00032702028 

26.1 

00 

43 

10.04 

+41 

16 

54.2 

2550 

9.6 

49 

00032702029 

27.2 

00 

43 

10.52 

+41 

17 

33.9 

2524 

10.2 

59 

00032702030 

28.4 

00 

43 

06.50 

+41 

17 

25.9 

2331 

10.2 

41 


In lBarnard et al.l (120061) we modeled the 2004 July 17 
XMM-Newton pn spectrum of XB158 with a blackbody 
and a power law, finding kT = 2.0±0.2 keV, the pho¬ 
ton index T = 2.0±0.3, \ 2 /dof = 19/19 and the 0.3-10 
keV luminosity was 1.5±0.6xl0 38 erg s -1 ; fitting a sin¬ 
gle power law emission model yielded a photon index of 
0.57±0.09, which is harder than any black hole spectrum 
(iRemillard & McClintocl3l2006l l. meaning that the accre- 
tor is a neutron star. We conducted three dimensional 
smoothed particle hydrodynamical (SPH) modeling, as¬ 
suming a 1.4 Mq NS, a 2.78 hr period, a total system 
mass of 1.8 Mq, and a luminosity ~80% of the Edding¬ 
ton limit. As a result, we estimated the disk precession 
peri od to be 29±1 t i mes t he orbital period, i.e. 81±3 hr. 

In IBarnard et al.1 (12012T) we examined Chandra obser¬ 
vations of ~30 M31 globular cluster X-ray sources span¬ 
ning ~12 years (89 ACIS and 45 HRC), including XB158. 
We observed 0.3-10 keV luminosities ~4-20xl0 37 erg 
s , and proposed that this is due to varying accretion 
rates over the disk precession cycle. 

Since we expected a super-orbital period for XB158 
over time-scales of a few days, we obtained a series of 30 
Swift observations, each with 2.5 ks exposure and spaced 
~1 day apart. In this work we present our analysis of 
these Swift results, and find evidence for a super-orbital 
period of days. 

2. OBSERVATIONS AND DATA ANALYSIS 

We made 30 Swift observations over 2013 February 
8 - March 9 (ObsIDs 0032702001-0032702030, PI R. 
Barnard). XB158 was one of two main targets for this 
survey, and the pointing was chosen to optimize the re¬ 
sults from both targets. However, the actual pointings 


varied significantly over the 30 observations. A journal of 
Swift observations is provided in Table Q] for each obser¬ 
vation we give the time relative to the first observation, 
pointing , XRT exposure, and off-axis angle, along with 
the net source counts. 

For each observation, we placed circular regions around 
XB158 and a suitable background region. We obtained 
the net source counts using the “Counts in regions” tool 
in the DS9 image viewer, and estimated the intensity by 
dividing the net counts by the exposure time. These data 
were obtained in order to determine the extent to which 
the varying off-axis angles affected our results. 

We also created spectra from the same extraction re¬ 
gions, created appropriate ancillary response files using 
XRTMKARF, and found the appropriate response file 
using the QUZCIF tool. None of the spectra were suit¬ 
able for free spectral fitting, hence we obtained lumi¬ 
nosity estimates by assuming a model obtained from 
previous observations. For Chandra ACIS observations 
with >200 net source photons, the mean line-of-siglit 
absorption (Ny) was 9±4xl0 20 atom cm -2 (x 2 /dof = 
7/11), and the mean power law index (T) was 0.52±0.03 
(y 2 /dof = 5/11). This is consistent with the best ab¬ 
sorbed power law fit to the 2004 July 17 XMM-Newton 
observation (Ah = 0.1, T = 0.57±0.09, x 2 /dof = 30/24 
IBarnard et al.l [20061 1. As we noted earlier, that XMM- 
Newton spectrum was best described by a blackbody + 
power law model, but neither the Chandra nor Swift 
spectra were sufficient to constrain the two-component 
emission model. 

However, we were able to estimate the luminosity for 
each observation by assuming a fixed emission model, al¬ 
lowing only the normalization to vary. We fitted each 
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XB158 Swift XRT 



Figure 1. 0.3-10 keV Swift luminosity light curve of XB158 for 

our 30 days of observations. 


spectrum using XSPEC 12.8.2b, fixing TVh = 9xl0 20 
atom cm -2 and T = 0.52, to find the normalization re¬ 
quired to make absorbed model intensity 1.00 count s — . 
We then calculated the unabsorbed flux for this model, 
allowing us to convert from intensity to flux. Multiplying 
the conversion by the background-subtracted intensity 
provided by XSPEC yielded the instrument-corrected, 
background-subtracted source flux. The luminosity was 
calcul ated f rom the flux assuming a distance of 780 kpc 
(iStanek fc Garnavichlll998 ). 

We fitted the lightcurve with constant and sinu¬ 
soidal components using the QDP program provided in 
HE ATO OLS, performing a simple search for periodic¬ 
ity. iScarglel (119821) created a periodogram suitable for 
unbinned data with a mean of zero that produces ex¬ 
actly equivalent results to such least-squares fitting of 
sinewaves, but also allows comparison of the best period 
with other periods. The likelihood that any peak in the 
periodogram is real is given by the false alarm probabil¬ 
ity (P), where a low value of P indicates that the peak 
is likely to be significant. If the highest frequency to 
be probed is N times higher than the lowest frequency, 
then the power, z , required for a false alarm probabil¬ 


ity P is given by z = — In 


1 - (1 ~P) 


1/N 


z ~ ln(AT/P) (Scarglc 1982). 


for small P, 


3. RESULTS 

We present our 30 day, 0.3-10 keV Swift XRT 
lightcurve of XB158 in Figure [Q We see that the 0.3- 
10 keV luminosity varied by a factor ~5; the luminos¬ 
ity dropped from ~2.3xl0 38 erg s -1 to ~5xl0 37 erg 
s -1 in ~2 days, assuming the mean Chandra absorbed 
power law model. We note that these Swift observa¬ 
tions are non-contiguous, spacing the 2.5 ks observing 
time over several hours; the low intensities (^0.005-0.025 
count s -1 ) meant that there was no appreciable variabil¬ 
ity within each observation. 

We find no evidence for a dependence of luminosity on 
off-axis angle; the luminosities for the observations when 
XB158 have the largest and smallest off-axis angles have 
consistent values, while observations at an off-axis an¬ 
gle of 8.2' resulted in a factor ~5 range in luminosity. 


Swift XRT lightcurve folded on 5.65 day period 



Figure 2. Folded 0.3-10 keV Swift XRT luminosity lightcurve, 
assuming a 5.65 day super-orbital period. 

The conversion factor for translating 1 count s -1 into 
0.3-10 keV flux ranged over 1.09-1.29x 10 -1 ° erg cm -2 
count -1 for all observations except the second one, where 
it was 1.79xlO -10 erg cm -2 count -1 . This ~10% varia¬ 
tion about the mean in instrumental correction is clearly 
not sufficient to account for the factor ~5 variation in 
luminosity. 

Fitting our lightcurve with a constant intensity yielded 
a best fit luminosity of ~1.3xl0 38 erg s -1 , with y 2 = 181 
for 29 degrees of freedom (dof). However, adding a si¬ 
nusoidal variation component yielded a much improved 
fit, with y 2 /dof = 43/26; for this model, the period is 
5.65±0.05 days, with an amplitude of 7.1±0.6xl0 37 erg 
s -1 around a mean luminosity of 1.43±0.04x 10 38 erg 
s -1 , and a phase of 88.1±0.7 degrees. All uncertainties 
in this work are quoted at the ler level. 

This sinusoidal variability yielded Ay 2 = 138 for 30 
bins, with 3 extra free parameters; F-testing showed that 
the probability for this improvement being due to chance 
was 3xl0 -8 , equivalent to a >5cr detection. In Figure [2] 
we show our Swift lightcurve folded on a 5.65 day period 
and fitted with the best fit sinusoid. 

We present our Lomb-Scargle periodogram for the 0.3- 
10 keV Swift XRT luminosity lightcurve of XB158 in Fig¬ 
ure [3l and also indicate the power required for false alarm 
probabilities P = 0.5, 0.05, and 0.005 for reference. We 
tested 30 frequencies, and oversampled each frequency 
by a factor of 50; this resulted in a periodogram cover¬ 
ing a wider range of periods than is interesting (up to 
~1500 days), so we show only part of the periodogram 
here. The periodogram shows a single strong peak, with 
the maximum power of 10.4 at a period of 5.60 days, 
corresponding to P = 0.0009; the range of periods which 
are detected at a >3cr level is 5.4-5.8 days. While there 
is a small peak at the 3 day period, P = 1 for this peak. 

4. DISCUSSION AND CONCLUSIONS 

XB158 is a high inclination X-ray binary associated 
with the M31 globular cluster B158. It exhibits deep 
intensity dips on a 2.78 hr period i n som e obse rvations 
but not others, prompting iBarnard et all (|2006H to sug¬ 
gest that the disk is precessing, caused by the “super¬ 
humping” phenomenon observed in low mass ratio sys¬ 
tems where the disk crosses the 3:1 resonance with the 
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Swift XRT LS periodogram for XB158 (P = 0.0009) 



Figure 3. Lomb-Scargle periodogram for our 0.3-10 keV Swift 
XRT luminosity lightcurve, identifying the power required for false 
alarm probabilites (P) = 0.5, 0.05, and 0.005. We see a peak at 
5.6 days with P = 0.0009, while 5.4-5.8 day periods are significant 
at a >3cr level. 

do nor star. 

iBarnard et al.l ( 2006 ) predicted a disk precession pe¬ 
riod of 29±1 times the orbital period, i.e. 81±3 hr, in¬ 
spiring a month of daily monitoring of the M31 central 
region by Swift. Fitting a sinusoid to the lightcurve re¬ 
vealed a 5.65±0.05 day super-orbital period (ler uncer¬ 
tainties); the 0.3-10 keV luminosity varied by a factor 
~5, which is consistent with the range of luminosities ob¬ 
served in the ACIS obs ervations o f our 13+ year Chandra 
monitoring campaign (Barnard et al.l 1 201 2f l. 

The peak of the Lomb-Scargle periodogram is at 5.6 
days, consistent with that obtained from least squares 
fitting of the lightcurve with a sinusoid. The suggested 
super-orbital period is ~ 70% longer than pre dicted by 
our 3D SPH simulations (IBarnard et al.1 12006 ). None of 
the authors of the current paper are experts in SPH; 
however, J. R. Murray stated in a private communica¬ 
tion that the longer than expected super-orbital period 
is likely due to the mass ratio of the donor to the accre- 
tor being lower than assumed. The mass of a Roche lobe 
filling main sequence star, m, may be approximated as 
m ~ 0.11 Ph r , where Py, ,. is the orbital period in hours 
(|Frank. King, fe Rainel120021 ). and we originally assumed 
a donor mass of ~0.30 M & . For the accretor, we as¬ 
sumed a 1.4 Mq neutron star. A power law emission 
model fit to the 2004 July 16 XMM-Newt on spectrum of 
B158 yielded a photon index of 0.57+0.09 (Barnard et al. 
2006), which is considerably harder than any spectrum 
emitted by a black hole binary. H owever, some n e utro n 
stars have masses >2M 0 (see e.g. iDemorest et ahll201(H ; 
iLvnch et al.ll2013h . and it is possible that XB158 contains 
a particularly massive neutron star. 


We noted in lBarnard et al i (i2012f ) that other XBs such 
as XB146 exhibited strong luminosity fluctuations be¬ 
tween fairly consistent maxima and minima during our 
Chandra monitoring observations, and suggested that 
this long-term behavior could be indicative of a short 
period / low mass ratio system. Our new findings sup¬ 
port this hypothesis. 
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